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Abstract 


Operations research methods are among the modern scientific methods that have occupied a prominent place 
among the mathematical methods used in planning and managing various economic and military activities. They 
have been able to help specialists in developing ideal plans in terms of costs, production, storage, or investment of 
human energies. One of its most important methods is the method Linear programming, which was built based on 
the sets of linear equations that represent the constraints for any linear model. Based on the methods for solving 
the systems of linear equations, researchers were able to prepare algorithms for solving linear models, such as the 
direct Simplex algorithm and its modifications. After the emergence of neutrosophic science, we found that 
research methods had to be reformulated. Operations using the concepts of this science, and as a basis and 
foundation for neutrosophic linear programming. In this research, we will reformulate the systems of linear 
equations and some methods for solving them using the concepts of neutrosophic to be a basis for any study 
presented in the field of neutrosophic linear programming. 


Keywords: Operations research; linear programming; systems of linear equations; neutrosophic science; systems 
of neutrosophic linear equations; methods for solving neutrosophic linear equations. 


1. Introduction: 

Through our endeavor to provide everything new and keep pace with scientific development and given the great 
importance of the linear programming method as one of the methods of operations research, we found it necessary 
to reformulate the systems of linear equations and some methods for solving them using the concepts of 
neutrosophic, the science in which research and studies that use its concepts topped the research list in most 
publishing houses. The reason for this is the accuracy of the results obtained in studies that use the concepts of this 
science, some research [1-13]. 


Operations research methods are among the modern scientific methods that have occupied a prominent place 
among the mathematical methods used in planning and managing various economic and military activities. They 
have been able to help specialists in developing ideal plans in terms of costs, production, storage, or investment of 
human energies. One of its most important methods is the method Linear programming, which relies on 
mathematical models to express any practical system from real life or a proposed idea. Mathematical models 
consist of an objective function and constraints. Classical operations research methods provided many algorithms 
that help to find the optimal solution for these models. The construction of these algorithms was based on sentences 
Linear equations. Based on the methods used to solve these sentences, such as the Jardin-Diver method and the 
Simplex method, the direct Simplex algorithm and its modifications were developed to solve linear models, and 
in view of the importance of the linear programming method and to keep pace with everything new in the field of 
scientific research and after the emergence of neutrosophic science, the science that caused a great revolution. In 
all fields of science, we found that it is necessary to reformulate the methods of operations research using the 
concepts of this science, the science that takes into account all the changes that may occur to the issue under study 
through the indeterminacy of the issue’s data. Therefore, we present in this research a study of sentences of linear 
equations and some methods for solving them using the concepts of neutrosophic science to be a basis for any 
study presented in the field of neutrosophic linear programming .we will take the matrix of constants on the right 
side of these equations as neutrosophic values. Of the form Nb; = b; + 6; where 6; is indeterminacy and can take 


one of the formsd; € {ui1, Hiz} or 6; € [Wir, Hi2] . , also aj; on the right side we will take it Neutrosophic values 
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ie., Nay; = aj; + yi; where y;jis indeterminacy and can take In one of the forms y;; € [Gijr Giza], or Vij € 
{9; jv Pi ja}; which helps us obtain more accurate solution results that take into account all conditions. 


2. Discussion: 
Based on the study mentioned in the classical references on the systems of linear equations and methods for 
solving them in some references [14-17]. 

We present the following study according to the concepts of neutrosophic science: 
The systems of linear equations in which the number of equations equals m and the number of variables in them 
equals n are given according to classical logic in the following general form: 


44X41 + Q42X2 + vom + AynXn = by 
24X41 + A22X2 Sie A2nXn = bo 
Am1X1 + Am2X2 t+ + AmnXn = Dm 
It is written in the following matrix form: 


A.X=B 
Where: 
ay Aq2 + Ain b, X 
ar Az? «Aon b X2 
A= B=| ? X= 
ami am? «Amn Din Xy 


Where aij and b; are real numbers for all values of i = 1,2,...,m andj = 1,2,...,n 
We distinguished three cases of the systems of linear equations: 
The first case: The number of equations is equal to the number of variables, i.c.,.m =n 
The second case: The number of equations is greater than the number of variables, i.c..m >n 
Third case: The number of equations is less than the number of variables, i.e., m <n 
Below we will present the systems of linear equations using the concepts of neutrosophic science, where we will 
take the real numbers aj and b; as neutrosophic numbers, that is, of the form Nb; and Naj,;, indefinite values. 
Perfectly determined, they can be any neighborhood of the real numbers a;; and b;, written in one of the forms. 
next: 

Na, = a,j + €; and Nb; = b; + uw; where &;; € [Asij Avis] or &; € {Arij Avis} or otherwise, then the systems of 
neutrosophic linear equations is written in the following form: 
The systems of neutrosophic linear equations in which the number of equations equals m and the number of 
variables in them equals nare given in the following general form: 


NQ44X1 + NQy2X2 +++ + NAyyX_ = Nb, 
NQy4X1 + NQg2X2 +++ + NAgyX, = Nb 


NamnX%1 + NAmnX2 +°°° + NAmnXm = Nb 


In the following matrix form: 


NA.X = NB 
Where: 
Na,, Na,..Nain Nb, Xy 
NA = Naz, Naz  NQon NB = Nb, X= X2 
Nain, NOgup N Genm Nbm Xn 


To find the general solution to the previous systems of equations, we examine them according to the three 
aforementioned cases: 


The first case: The number of equations is equal to the number of variables, i.e., m =n. 
We write the systems of equations as follows: 

NQ44X1 + NQy2X2 +++ + NQyyX_ = Nb, 

NQy4X1 + NQg2X2 +++ + NAgyX, = Nb 
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NQy1X1 + NAngX2 + 2+ + NAgnXn = Nby 


In matrix form: 


Where: 


NA.X =NB 
Na,, Naj,..Nain Nb, xy 
NA _ Naz, Naz2 Non NB = Nb, xX a X2 
NQn1  Nanz »NQnn Nb, Xn 


The proverbs matrix is a square matrix whose determinant is Ay= |NA| 
Here we distinguish two cases: 


1- 
a. 


Ay= 0 . This case results in two cases: 
If Ay= 0 and Ay, # 0 where Ay, is the determinant resulting from the determinant of the matrix of 
Jj ij 


proverbs Ay after replacing the column containing the proverbs of the unknown x; with the column of 
constants (the values on the side The second of the equations) then the systems have no solution. 


If Ay= 0 and Ay, = 0, this means that the systems of equations are not linearly independent, meaning 
ij 


that some of them are linearly related to each other. To address this case, we delete one of the two linearly 
related equations, thus the number of equations decreases and becomes m’where m’ = m — 1 andm’ < 
n, which is identical to the second case that will be dealt with later. 
If Ay# 0, that is, the systems of equations are linearly independent and the systems have a single solution, 
which can be found in several ways. In this research, we study Gaussian- Jordan method, which is the 
basis for the direct simplex algorithm that we use to obtain the optimal solution for linear models. 
Gaussian- Jordan method for solving systems of neutrosophic linear equations in which m = n: 
To clarify the mathematical basis of this method, we write the equations in the following matrix form: 
Nay, Na, ..Nayy xy Nb, 
Naz, Naz ...NQ2n iG _ Nb, (1) 
Nay, = Nan2+»NQnn 1 Xn Nby, 
Or in the following abbreviated form: 
NA.X=NB (2) 
Since Ay= |NA| ¥ 0, this means that the matrix NA has an inverse that is NA~*. We multiply both 
sides of equation (2) by NA~*and we find: 
NA~*.(NA.X) = NA7*.NB 
Hence, we get: 
I.X=NB' 
Which is written in the following detailed form: 
10 0..0] [% Nb, 
0 1 sca is _ |Nb, (3) 
0 0 0..1]1 LXn Nby, 
This process is the basis of the Gaussian- Jordan method for solving a system of linear equations. In 
order to convert Figure (1) to Figure (2), we follow the following steps: 
1- We express Figure (1) in the following table: 
Table No. (1): Table of equations 


Variables x4 Xo oa ae NB 
Equations 

1 Nay, Nay oT Nain Nb, 

2 Nay, Nag Be Nazn Nb, 

n Nani Nan2 ies Nann Nby 


2- We convert the matrix NA to the unit matrix J by processing the lines of the table so that we make 
all non-diagonal elements in all its rows equal to zero and the diagonal elements equal to one. When 
we want to remove the variable x, from the equation t, we follow these steps: 

a- We divide all the elements of row t in which we want to make x, equal to one by Na,,, so x, becomes 
equal to one and the other expressions change. 

b- We set all elements of the column with x, (except row t) equal to zero. 

c- We calculate the rest of the elements of the new table from the following two relationships: 
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aati) _ NajjNa,, — Na;,Naz; 

7 Nats 

Nb;Na,, — Naj;Nb;, (4) 
Nats 


The element Na,, is called the pivot element 
From the previous processing we get the following table: 
Table No. (2) Final solution table 


Variables | Nx, Nx eed NXn NB' 
Equations 
1 1 0 ees 0 Nb; 
2 0 1 = 0 Nb; 
n 0 0 ai 1 Nby 
Through the table, the systems of linear equations are written in the following matrix form: 
I.NX = NB' 
10 0..0] [% Nb, 
0 1 0..0] |*2 Nb; 


= => 


00 0..1) Lem] Lb, 


Nx, Nb, 
Nx2 = Nb3 
NXn Nb, 


=> Nx, = Nb, ,Nx, = Nbj,...,Nx, = Nb, 
The second case: 
The number of equations is greater than the number of variables, i.e., m > n- 


In this case, we form a new system from the set of equations in which the number of equations is equal to the 
number of variables by excluding a number of equations of m — n. Then we solve the new systems as we did in 
the first case and replace the resulting solution in the equations that were excluded to ensure that they are satisfied. 


The third case: The number of variables is greater than the number of equations, i.e., m < n- 
In this case, we are faced with a set of equations of the following form: 
NQ44X1 + NQy2X2q +21 + NA Xm +211 + NA yX_ = Nb, 
NQg14X1 + NQg9Xq ++ + NQgmXm +11 + NAgnX_ = Nb 


eee ee eee eee ee eee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee (5) 
NamixX1 + Nam2X2 + °° + NQmmXm ++: + NAmnXn = Nby 
Which is written in the following matrix form: 

x41 
x2 

Nay, = Nayz . NQqyy we vs Nain Nb, 
NA= Nay, NQg2 .» NQgyy we Nain NB = Nb, x =|X%m 

NQm1  NQyn2 +» NQmm - Namn Nb, 
Xn 


In the following brief form: 
NAmmay -Xavay = NBanay (6) 
1- We partition the matrix NA (mn) into two matrices: 
a- A square matrix of rank (m.m) and we denote it as NC(m m)- 
b- And a rectangular matrix of rank (m.n — m) and we denote it NDann—m) 
2- We partition the column matrix X( 1) into two matrices X' (m1) and X : (n—m.1)- 
Then the systems of equations (5) are written in the following matrix form: 


= NBom.1) (7) 


NConm): X' (m1) ot NDann-m)- x’ (n—m.1) = NBanay 
From them we find that: 
NConm): ent) = NBs) ~ NDomn-—m): x (n-—m.1) (8) 
Assuming that |NC| # 0, we multiply both sides in relation (8) by NC~? and we find: 
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NC71.NC.X' = NC71.(NB — ND.X’) 
I.X' = NC-1.NB — NC~1.ND.X") (9) 
Assuming that NC~*.NB = NB' and NC~*.ND = NDonn—myWe find that: 


10 0O...0 xy Nb, Ndi, Nd, -Ndi@—m) Xm41 
0 1 0..0) |}%2]_|NbZ]_ | Ndi, Ndjy «Nim | |*m+2 (10) 
0 0 0..11 Xm Nby Ndny = Ndi -- Ndinmn—m) Xn 


It is transformed into a set of linear equations as follows: 
Nx, = Nb, — (NdyyXm41 + NdipXma2 + + Ndia—myXn) 
Nx2 = Nbz — (Nd31%m41 + Ndz2Xmy2 + °° + Nd3(n—m)Xn) 


NXm = Nbn aad (NdiniXm+41 + Ndin2X%m+2 ae tes Ndin(n—m)*n) 
This means that we were able to calculate m in terms of (n — m)in terms of (n — M), Xm41,%Xm4+2) ++» Xn, We note 
that the values of the variables x1, X2,...,Xm, it relates to the values taken by the variables Xm44,Xm42)++)Xn, OF 
in other words, what we give to the variables X44, Xm42) ++», and that for every sentence of Values such as 
Bin+1» Pm-+2) ++» Py for these variables we get a set of values for the variables x, X3, ..., Xm iS: 
Nx, = Nb, = (Ndi1Bm41 + Nd12Bm+2 ap eet Ndi (n-m)Pn) 
Nx, = Nb, aa (Nd31Bm41 + Ndj2Bm+2 toe Nd3(n-m)Pn) 
NXm = Nb = (NdintBm+1 +e Ndin2Bm+2 ab rtek Ndinci-myPn) 
Thus, we obtain a solution that includes all the variables of sentence (5) 
The solution is arranged as follows: 
(B1, B2, sey Bm Bm+v Bm+2) eae Bn) 
But since the variables x41,%Xm+42) +, %,can take an infinite number of qualitative values (even if they are 
restricted by certain conditions), we obtain an infinite number of corresponding values for the variables 
Ky Kop snp Ais 
Therefore, if |NC| # 0, then the set of equations (5) has an infinite number of acceptable solutions of the form: 
(Xp) X20 0 Xm X41 Xme4a0 Xn) 
Thus, we obtain a solution that includes all variables of the sentence, which is the ordered solution: 
(Pr, B2, ey Bm Bm+v Bm+2) ee Bn) 
3. Basic solutions 
Since sentence (5) has an infinite number of acceptable solutions, we will try to limit ourselves to a limited number 
of them by setting the variables X41,%Xm+42,+.,Xp, equal to zero. Then sentence (9) takes the following form: 
10 O...0 xX Nb; 
xX y 
0 1 0..0 Bas Nb, (11) 
0 0 0..11 L%m Nby 
From it we get: 


X, =Nb, ,X2 =Nb3,..,Xm = Ndi, 
So, the complete solution is: 
(Nb;, Nb3, ..., Nbj,, 0,0, ...,0) 
We call this solution the basic solution because it is attributed to the rule with single normal vectors in the space 
R™ as follows: 
1 0 0 
0 1 0 
0 0 0 
0 0 0 
0 0 m 
The set of vectors €,,@2,..-,€m form a rule because they are linearly independent, and the vector NB’ can be 
expressed in terms of it using the factorials x1, x2, ..., Xj, as follows: 
NB’ = €4X1 + 9X2 ++ + em@Xm 
We call the variablesx,, Xz, ...,%m, basic variables and we call other variables Xj744, Xm42) ++ )Xn 
free or non-basic variables because they take qualitative values. 
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The process of choosing the variables x1, Xz, ..., X,;,, to be basic variables is a random process, as we can form other 


basic solutions, knowing that the possibilities available to obtain basic solutions are: 
n! 
m — 
on = m! (n —m)! 
It is a finite number of infinite acceptable solutions. 
Example 1: 
Find the joint solution of the following two linear equations: 

2x, + 7X2 + 3x3 4+ 2x, = [2,5] 

3x, + 9x, + 4x34 x, = [3,7] 

X, + 5x, + 3x, +4x, = [4,8] 
In the set of equations, the number of variables is n = 4 and the number of equations is m = 3. Therefore, the 
number of basic variables is equal to 3 and the number of non-basic free variables is n — m = 1. The number of 
possible solutions is calculated from the relationship: 


n! 
m _ 
Cn = m!(n —m)! 


4) 
Ci = aaa) 
Write as follows: 
(X1,%X2,%3,0) , (%1,%2,0, x4) , (%1,0,%3, X4), (0, X2, X3, X4) 
To obtain these solutions, we write the systems of equations in the following form: 
2x1, + 7X2 + 3x3 = [2,5] — 2x, 
3x, + 9x2 + 4x3 4X4 = [3,7] —X4 

X, + 5x2 + 3x3 = [4,8] — 4x, 

The previous sentence is written in the following matrix form: 


2 9 a F [2,5] 
F 9 4. = |[3,7]} - -[i]. 04 [x4] (*) 
1 5 3 fe [4,8] 
2 7 3 [2,5] 
c=, 9 4 xX'= wy NB = |[3,7] D= ; x= (ea ] 
1 5 3 x3 [4,8] 4 
We calculate the determinant |C| . 
We find: 
2 7 3 
|C| = 9 44=-340 
1 5 3 
2 7 3 
To find the solutions, we find the reciprocal of the matrix, C =|3 9 ] 
15 3 
We find: 
—7 —1 
= =o 
C+=]5 —1 
3 - 3 
—2 1 1 
We compensate in the relationship: 
NC71.NC.X' = NC~1.(NB — ND.X’) 
We get: 
—7 —1 —7 = 
3 * 3le73 ls 2 SI /f25) p 
5 -1 iF 9 ‘|= 5 —1].| | [3,7] - |i} a 
3 ~ gslh 5 3) |3 7 [4,8] 
—2 1 1 —2 
—7 —1 —7 
10 o1fNay [3 % 3B) [25] |= 
lo 1 ollvs]—[5 —1}.][3,7]]-| 5 le 
oo allwx,) | 3 -! ZF] I [48] 3 = 
—2 1 1 - 
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It is transformed into the following systems of equations: 


lo >] 16 
2° ¢ 
Nx, =| _ [1.5] — . [x4] 
s 1 
[3,5] 
Setting the free variable x, equal to zero we get: 
—1 
0,— 
5 
Nx, = 4 
a | ee 
[5 
[3,5] 


—1 4 
xy = lo, x. = = Fee = [3,5] 


Thus, we obtain the first neutrosophic basic solution, which is: 


(x1, Xz, X3,0) = (fo :— [1.5] , [3,5], 0) 


In the same way we obtain other basic solutions. 


4. Dissolved basic solutions: 
The base solution is a degenerate and invalid solution if in the final result we obtain a value of zero for the variables 
that we chose as the base. 
Gaussian- Jordan method for solving a set of linear equations in which m < n: 
Based on the previous mathematical principles, the basic steps of the Gaussian- Jordan method are as follows: 
1- We write the systems of equations (5) in the following matrix form: 
1.X'+NC71*.D.X° = NC1.NB = NB' 


(1, NC-.ND]. [| = NB’ (12) 


Which is written in the following detailed form: 
10 0..0 Ndi, Nd}, ..Ndi cm) X1 Nb} 
0 1. 0..0 Ndi, Ndjz ...Ndim—m) X2| _ | Nb, (13) 
0 0 0..1 Nding Ndinz « Ndingn—my} Vn Nby, 
The transition from Figure (5) to Figure (12) is done by following the same steps that we mentioned in the previous 
paragraph, but this method does not give us a basic solution unless we set all the free variables equal to zero. If we 
do that, we only get the first solution, and to get all the solutions, we do the following steps: 
a. We organize the following table: 
Table No. (3) The first table for the Gaussian- Jordan method 


Variables x4 X2 ee Kon. Xin Xm+2 bee Xn NB 
Equations 
ee 
1 Q41 a42 tee aim Qim+1 Aim+2 oo Gin Nb, 
t 
2 a21 a22 tt a2m Agm+1 Agm+2 ve Gon Nb, 
y 
m Am1 Am2 nes Amm. Amm+1 Amm+2 ae Amn Nbn 


b. We find the identity matrix I,,.. by processing the rows of the previous table in the same way that was 
explained in the previous paragraph. This is done by specifying the variables that will be entered into the 
base and let them be x1, X2,...,X,. AS a result of this processing, we get the following table: 


Table No. (4): Table of the first basic solution 


Variables Xx X2 aoe | ad Xoned Xm+2 a Lon NB' 
Equations 

1 1 Ol eee Ndi, | Ndi, |... | Ndin om Nbt 

2 0 7 la ow Nd, | Ndb. |. | Ndbn—m Nb! 

m 0 o |... | 1 | Ndi, | Ndb, ae. oo Nbi, 
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c. Setting all the free variables in Table (4) equal to zero, we obtain the following first basic solution: 
(Nb;, Nb3, ..., Nb}, 0,0, ...,0) 

d. To obtain a second basic solution, we replace one of the basic variables, say x,,, with one of the non- 
basic variables x41, by selecting the appropriate pivot element, and here it is Nq,,,. We work to delete 
Xm-+1 from all equations except equation m .We set the coefficient of this variable in this equation equal 
to one. We perform the appropriate calculations through the two relations (4). We solve the following 
second basic solution: 

(Nbi, Nb5, ..., Nby_1,0, NDin4 4, 0, 0) 

To obtain other basic solutions, we repeat what was stated in step (d). 

Non-negative basic solutions: 

If all or some of the variables are conditioned to be non-negative, then some of the basic solutions are unacceptable 

because they violate the condition. In such a case, we have to look for positive basic solutions from among the 

basic solutions. 

Due to the difficulty of applying the method mentioned in the example, especially in the case that includes a large 

number of variables, the Gaussian- Jordan method was developed so that positive solutions are directly obtained. 

The new method was called the simplex method, which is carried out according to the following steps: 

The simplex method for finding non-negative basic solutions to a system of linear equations in which m < 

n: 

In the systems of equations (5): 

1- Wemake all elements of the constant’s column NB on the second side of the equations non-negative, 
by multiplying the equation whose second side is negative by (-1) 
2- We put the coefficients of the new systems in a table. 
3- We form a rule consisting of m variables by choosing the variable that we want to enter into the rule, 
for example, x,, and then we calculate the index. 
ga - _ Nb, 
oe Najs] Nats 

We call the element Na;,, the pivot element, we delete the variable x, from all equations according 

to the Gaussian- Jordan method, except for the equation t, in which its coefficient is equal to one. 

We repeat the previous step until we form a base consisting of m variables. 

4- Setting the non-basic variables equal to zero we obtain the following non-negative basic solution: 

(Nb,, Nb3, ..., NBy_1, 0, Nbin 44, 0, ...,0) 

5- To obtain other new non-negative basic solutions, we choose one of the variables to be a basic 
variable, then we determine the pivot element and repeat the work we did for the previously 
mentioned variable x,. We obtain a new non-negative basic solution, and so we continue working 
until we obtain all non-negative basic solutions. 


> 0; Na;, > 0, Nb; > 0 


We explain the above through the following example: 
Example 2: 
xX, — 3x, + 2x, = —[1,3] 
Xp + 2x, — 3x5 = [2,8] 


We multiply the first equation by (-1) until the condition Nb; > 0 is met, we obtain the following new systems: 
—x, — 3x3 — 2x, = [1,3] 
X2 + 2x, — 3x5 = [2,8] 
The stopping criterion is when we do not find a free column that we did not use for switching that contains a 
positive element (at least one), meaning that all the elements of the free columns that were not used during the 
swap are negative values. 
In the systems of equations, the number of variables is n = 5 and the number of equations is 
m = 2 Therefore, the number of basic variables is equal to 2 and the number of non-basic free variables is n — 
m = 3. The number of possible solutions is calculated from the relationship: 


n! 
cm = —____ 
”" — m!(n-—m)! 
1.€., 
; 5! 
Cs = ———— = 10 
5 21(5 — 2)! 


Write as follows: 
(X41, Xz, 0,0,0) , (x, 0, x3, 0,0), (x3, 0,0, x4, 0), (x4, 0,0,0, x. ), (0, x2, x3, 0,0) 
(0, x2,0, x4, 0), (0, x2, 0,0, xs), (0,0, x3, x4, 0), (0,0, x3, 0, x5), (0,0,0, x4, x5) 
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To obtain these solutions, we organize the following table: 


Table No. (5): The first table for the simplex method 


Variables x4 X2 X3 X4 Xs NB 
Equations 

1 -1 0 0 3 =2 [1,3] 

2 0 1 0 2 —3 [2,8] 


To find a basic solution to the set of equations, we choose a variable, for example x4, to be a basic variable, and 
to determine the appropriate anchor element, we calculate the index: 


i 3; 2 3 


us 


p= tn] = en 4, a] = 


That is, the fulcrum is a,, = 3. By performing the necessary calculations to delete the variable x, from the two 
equations, we obtain the following table: 


Table No. (6) The second table for the simplex method 


Variables x4 X2 X3 X4 Xs NB' 
Equations 
X4 1 0 0 1 2 1 P 
- =| Be 
2 2 iL 0 0 5 F 6 
3 3 3” 


We choose another variable to be a basic variable. We note that the variable x2 is ready to be a basic variable, 
and thus we get the following table: 


Table No. (7): Final solution table 


Variables X41 Xz X3 X4 Xs NB' 
Equations 
Xq 1 0 0 1 2 E . 
3 3 | Bl 
X2 2 1 0 0 5 E 6 
3 3 3’ 


Thus, we obtain a base consisting of the variables x, ,x,. We set the free variables equal to zero, and we obtain 
the following non-negative neutrosophic basic solution: 


(0.55 4].0.5-1)-0) 


To obtain other solutions, we repeat what we did to determine the previous solution. 


5. Conclusion and results 

As a basis for neutrosophic linear programming, we presented in this research a study of the sets of neutrosophic 
linear equations, and the Gaussian- Jordan method, which is considered the mathematical basis for the simplex 
method used to find positive basic solutions, which can be used when there are restrictions on some or all of the 
variables to be positive values, which in turn was the basis for the method. Direct simplex used to find the optimal 
solution for linear models. Through the examples that we presented on the systems of neutrosophic equations, we 
arrived at basic neutrosophic solutions that express unspecified values. Such sentences can be used in cases where 
the data provided to the systems that operate according to these systems of equations is subject to change. Here 
we can benefit from the margin of freedom offered by neutrosophic values. 
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